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Abstract 

The understanding of inactivation radiobiological mechanism in individual cells 
is important when from one side the application of ionizing radiation to tumour 
I ^ therapy and from the other side the protection against radiation are to be effectively 

improved. Important part of this mechanism is double-strand break (DSB) forma- 
tion; these DSB may be formed directly by impacting ionizing particles or indirectly 
by different secondary radicals. The latter kind of formation is much more frequent 
when cells contain normal water content. Mathematical model of the corresponding 
chemical stage will be presented with the aim to demonstrate how the individual 
radicals, but also other chemical agents present in the water during irradiation may 
influence DSB formation. 
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o | Introduction 

Ionizing radiation represents one of the main methods of tumour treatment. The 
aim is to destroy tumour focus without damaging irreversibly surrounding tissues. 
To fulfill this goal it is necessary to understand well the radiobiological processes 
running in individual living cells after the impact of a ionizing particle. The given 
results will be helpful also as to radiation protection. 

According to contemporary knowledge the radiobiological processes in a normal 
t-H cell (with usual water content) may be divided into three phases: physical (energy 

transfer from radiation to medium and the formation of water radicals), chemical 
(reactions initiated by radicals and formation of damages in chromosomal DNA) and 
biological (reaction of individual cells to the damage ending or by damage repair or 
by cell inactivation) . The first two phases are finished in a very short time interval 
(fraction of second) so that at dose-rates applied to in radiotherapy treatment (or 
playing role in radiation protection) the effects of individual ionizing particles occur 
practically quite independently In the radiotherapy the total dose is delivered in 
a rather short time so that the total DNA damage gathered during one irradiation 
may be regarded as the starting point of the biological phase ending after much 
longer time (at least ten minutes or longer). 

It is possible to say that the final result of biological phase depends on the 
number and distribution of individual DSB (double strand breaks) in the system of 
chromosomes. Thus the probability of DSB formation by one ionizing particle (or 
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photon) may be taken as the basic parameter determining inactivation effect, fn the 
following we will try to analyze the importance of individual processes running in 
the chemical phase with the help of a corresponding mathematical model. 

However, it is necessary to distinguish between the ionizing particles exhibiting 
different ionization densities. Photons, accelerated electrons, and in principle also 
protons belong to low-density or low-LET (linear energy transfer) radiation kinds 
when the cell being hit by one particle only are not inactivated as a rule. On the 
other side, the single accelerated heavier atom nuclei (high-LET radiation) inactivate 
cells having be hit by their Bragg peaks. As to the formation of DSB by low-LET 
radiation it may be rather strongly influenced by the presence of other (radiosensi- 
tive or radioprotective) chemical agents. And only rather sophisticated (complex) 
mathematical models may be helpful in understanding better corresponding mech- 
anism. 

Chemical phase and DSB formation 

We have already mentioned that the inactivation effect is approximately proportional 
to the efficiency of DSB formation. However, it depends also on the distribution of 
DSB in chromosomal system, which may be very different at a given dose for divers 
radiation types. The efficiency may be studied by measuring the survival curves, 
i.e. the ratio of surviving (not inactivated) cells for different radiation doses. And 
it is also possible to establish experimentally the numbers of DSB in the given 
chromosomal DNA under the same conditions. The latter method enables then 
to establish some important characteristics also outside the biological material for 
corresponding DNA molecules dissolved in water. 

One may assume that in DNA dissolved in water (i.e., also in normally living 
cells) practically all SSB and DSB are formed in indirect way by radicals formed in 
water. A direct collision of an individual ionizing particle with a given DNA may 
by practically neglected. The radical clusters of different sizes are formed along 
particle tracks. However, the most of them are so small that only one SSB may 
be formed with a certain probability. As to the low-LET radiation they are only 
clusters formed by track ends of electrons of greater original energy that may form a 
DSB with sufficient probability. Only heavier ionizing particles in their Bragg peaks 
are more efficient. The conclusion that DSB is practically always formed by a single 
radical cluster is strongly supported by the fact that the number of DSB increases 
linearly with increasing dose pQ. 

And we should ask what is the size of clusters that may form a DSB and how 
is the DSB formation influenced by chemical mechanism in individual clusters. Our 
goal is to study the reactions running in radical clusters and to analyze how their 
final effect depends on their size and might be influenced by the presence of other 
chemical agents; and what is the average size of radicals forming DSB for low- 
LET radiation. The problem may be hardly solved without the help of suitable 
mathematical models and computer simulations of the corresponding processes. 

After their impact the ionizing particles transfer their energy to the water medium 
and form clusters of water radicals (eventually, also of other ones) when other cor- 
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responding chemical agents are present, too. Oxygen radicals are usually involved 
under standard aerobic conditions. And it is known that the radiobiological effect 
increases with oxygen content, which should be accompanied by the corresponding 
increase of SSB and DSB numbers. It is important in radiotherapy as for low-LET 
radiation the effect in anaerobic tumor cells diminishes. 

Irradiating the cells with the normal content of water the chemical stage of 
radiobiological mechanism may have significant influence on the final effect. The 
following processes are then running in the corresponding radical clusters: 

(i) the free radicals having been formed by ionizing particle may mutually react 
and may recombine and form new chemical substances that react further; 

(ii) the radical clusters diffuse into surrounding medium, the density of radicals 
diminishes and their number gradually decreases; 

(iii) all substances present in the cluster (not only radicals) may take part in che- 
mical reactions and other agents (including radicals) may be formed; 

(iv) different agents (mainly chemically active radicals) present in the cluster may 
react with cell DNA and the biologically efficient DSB in chromosomal DNA 
may be formed. 

And it is evident that the efficiency of DSB formation will depend not only on the 
size of individual clusters, but also on all chemical agents present in them during 
irradiation, important role being played also by the oxygen in aerobic conditions. It 
is also the diffusion of individual clusters that may play an important role in final 
effect. And it is practically evident that rather complex mathematical models arc 
needed when we are to understand the efficiency of individual radicals in SSB (and 
DSB) formation. 

We have proposed a corresponding mathematical model involving the influence 
of chemical reactions as well as of diffusion; its preliminary version having been 
presented in Ref. [2] . We have demonstrated its application using the experimental 
data obtained by Block and Loman [3] in irradiating DNA molecules of bacteriofague 
<i>X174 dissolved in water by Co gamma radiation (of energy cca l,25MeV). The 
frequency of SSB (and DSB) formation in dependence on oxygen content in the 
mixture with N2 and N2O were determined in this experiment. 

Processes running in chemical phase 

The chemical processes running in water after the impact of ionizing particles has 
been described e.g. in the paper of Sauer and Schmidt [4]. Here a short survey of 
them will be given and their role in DSB formation will be mentioned. 

One may expect that they are mainly water radicals of H* , OH* and e~ q that 
form individual SSB, and also HO* radical if oxygen is present. And one may 
ask what is the average size of the radical clusters that would be consistent with 
frequencies of DSB established experimentally under different conditions. While 
HO* radical is formed in secondary reactions with oxygen all water radicals are 



3 



formed during the impact of ionizing particles. First of all it is the ionization of 
water molecule (ionization potential 12,56 ev): 

H 2 0^ H 2 0+ + e"; (1) 

H 2 + ion may decay and produce OH' radicals: 

H 2 + ^H + + OH'. (2) 

If the energy of a photon is approximately 7 eV the water molecule is excited and 
the following dissociation leads again to the formation of water radicals: 

H 2 -> H 2 0* -> H' + OH'. (3) 



Reactions 


Rate constants 

(dm^mor 1 ^- 1 ) 


1. 


H' + H* 




H 2 


1 ■ 10 1U 


2. 


Caq + H' 




H 2 + OH~ 


2.5 • 10 10 


3. 


&aq \ &aq 




H 2 + 20H~ 


6 ■ 10 10 


4. 


e- + OH' 




OH~ + H 2 


3 • 10 10 


5. 


H' + OH' 




H 2 


2.4 ■ 10 10 


6. 


OH' + OH' 




H2O2 


4 • 10 9 


7. 


H 3 0+ + e- 




H* + H2O 


2.3 • 10 10 


8. 


H 3 0+ + OH- 




H2O 


3 ■ 10 10 


9. 


H' + H2O2 




H2O + OH' 


1 • 10 8 


10. 


e~ q + H2O2 




OH' + OH~ 


1.2 ■ 10 10 


11. 


OH' + H2O2 




H 2 + HO' 2 


5- 10 7 


12. 


OH' + H 2 




H 2 + H' 


6- 10 7 


13. 


HO\ + H' 




H2O2 


1 ■ 10 10 


14. 


eaq + O2 




O- + H 2 


1.9- 10 10 


15. 


HOI + OH' 




H 2 + O2 


1 • 10 10 


16. 


HOI + HOI 




H2O2 + O2 


2 • 10 6 


17. 


H' +O2 




HOI 


1 ■ 10 10 


18. 


O2 + H 3 0+ 




H0' 2 


3 • 10 10 


19. 


H 2 




H 3 0+ + OH- 


5.5- 10" 6 


20. 


HO' 2 




H s O+ + 2 - 


1 • 10 6 



Table 1: Recombination reactions going in water after radiation impact and their 
rate constants [5] 

The electron released in reaction ([TJ may bind to a water molecule 

e~ + (H 2 0) n - e" (4) 

and electron e~ q contributes to production of hydrogen radicals: 

e~ q + H 2 O^OH- + H'. (5) 

Hydrated electrons e~ q may exist in this form relatively long; they may diffuse to 
greater distances and react with DNA molecules. 
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The effect of ionizing radiation is known to be modified if oxygen is present in the 
solution. It represents a source of some radicals which are assumed to be aggressive. 
They may arise by the reaction 



2 



o' + o- 



(6) 



minimum electron energy ~ 4 eV, maximum gain at 8 eV. These radicals react in 
water medium and new radicals are formed: 



O* + H 2 

O- + H 2 



H' + HO', 
H' + HO^. 



(7) 
(8) 



Oxygen radicals may be involved or formed further in the reactions 



2 + H 3 + 
H' + 2 
HO' 2 



2 , 
HO', 

HOI 

H.O+ 



(9) 

(10) 

(11) 
(12) 



At lower concentrations of 2 (at pH7) the dissociation (12 1 is preferred; at higher 



concentrations it is the formation of radicals HO' according to Eq. (10 1 which is 
preferred [5]. 

Thus, the actual content of radicals in an individual cluster depends also on oxy- 
gen concentration; the radicals H0 2 (resp. O') may be involved in DSB formation 
in addition to H' and OH'. Possible recombination processes and other reactions 
are summarized in Table [l] 



Mathematical modeling of the chemical mechanism 

As already mentioned the mathematical model proposed by us earlier (see [5]) has 
been applied to fitting the experimentally determined numbers of DSB for different 
concentrations of molecular oxygen in the mixtures of 2 with N 2 and N 2 in water 
solution of DNA. The proposed model was very simple. And we have succeeded 
in reproducing the data when the reaction rates (and also their mutual ratios) of 
involved chemical species differed rather significantly from those shown in Table 1. 

In the following we shall apply the model only to one set of experimental data 
(for the mixture of 2 with N 2 only), but we will broaden its basis. First of all, the 
model will be partially generalized and the number of recombination reactions will 
be extended. However, on the other side some limiting conditions will be imposed 
on the parameters that were left earlier quite free: 

- the diffusion coefficients for individual radicals will be taken from the literature 
and fixed (see Table 2); 

- it will be required for the ratios of the most reaction rates to correspond to 
experimentally established values (i.e., for mutual reactions of substances in clusters 
as well as for reaction of radicals with DNA) ; 
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Substance 


Diffusion coefficient 
(nm 2 .ns' 1 ^ 


Number 
of species 


Designation of 
dif. coefficients 


1. H' 


7.0 


N H 


D H 


2. OH' 


2.8 


Noh 


Doh 




4.5 


N e 


D e 


4. HOI 


2.0 


No 


Do 


5. H 3 0+ 


9.0 


N Ha o 


Dh 3 o 


6. 2 " 


2.1 


N o- 


D o- 



Tabic 2: Diffusion coefficients 



- and also for original number of radicals in corresponding cluster it will be 
required to correspond to efficiency of reactions responsible for their formation. 
Only several parameters (as will be mentioned later) will be left quite free. 

The mathematical model of corresponding chemical mechanism will start from 
the assumption of the existence of a hypothetical average cluster system character- 
ized by nonhomogeneous concentrations of individual species (see [7]). Macroscopic 
laws will be used to describe the diffusion of radiation-induced objects and the con- 
centration changes due to different chemical reactions. Such a situation may be 
described by the following set of coupled differential equations: 

^ = DiV 2 a ~ a £ kijCj + H k jk c i c k, (13) 

where Di are diffusion coefficients and Ci - spatially dependent concentrations of 
species i; kij are rate constant of reactions between species i and j. The first term 
of the right-hand side of the equation represent the diffusive contribution to the 
evolution of c% , while the second and third terms represent removals and productions 
of the i substances by chemical reactions. The set of such coupled equations may be 
solved cither by analytic approximation or numerically. 

We will consider the kinetics of low-LET radiolysis. The system will be described 
in terms of a spherically symmetric typical spur (cluster). We can then substitute 



equations (13 1 by equations, where diffusion processes are considered as spherically 
symmetrical 



j,k^i 



r denoting the distance from the cluster center. 

The most frequently used analytic models for the fast kinetics in radiolysis are 
based on approximation suggested by Jaffe (see [H]); they are known as prescribed 
diffusion. The initial spatial distributions of the radiation-induced particles, their 
concentration profiles, are assumed to be Gaussian, and the kinetic analysis invokes 
the approximation that reactions affect only the numbers of particles and not the 
form of the nonhomogeneous spatial profiles, which are therefore always Gaussian. 

One could solve the given system of partial differential equations numerically but 
it would be rather time demanding to determine the parameters of mathematical 
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model with the help of an optimization procedure (see e.g. [5]) to fit the used 
experimental data. Therefore, we have made use of the method which enables to 
substitute the solving of the system of partial differential equations by the solving 
of the system of ordinary differential equations, which is much less time consuming. 

The concentration profiles of species i due to diffusion process only may be ob- 
tained by solving Eq. ( 14) containing only the first term on the right-hand side; one 
obtains 



i(r,t) 



Ni 



-exp 



4Di* 



(15) 



where N% is the (initial) number of species i in the cluster. This concentration 
dependence corresponds to the Gaussian profile. Consequently, we can define the 
average concentration of species i as 



Ci(t) 



Ni 

Vi(ty 



(16) 



where Ni is the number of species i and Vi (t) is its time dependent average cluster 
volume at time t. The total number of species Ni does not change by influence of 



diffusion process. When we calculate the derivative of Eq. ( 16 1 we obtain 



dt 



Cj dVj 
Vi dt 



(17) 



which expresses the influence of diffusion process only. 



We can substitute now the first term of the right-hand side of the equation ( 14 1 



that expresses the influence of diffusion by the expession (17 1, which means that 



the diffusion of species i is expressed by corresponding time dependent volume Vi. 
Consequently, it is possible to substitute the system of partial differential equations 
by the system of ordinary differential equations 



dcj 
dt 



Cj dVj 
V dt 



X! k jk c 3 c k- 



(18) 



It means that in the fully described system also the numbers N decrease with 
time (see Eq. ( 16 1. And it is advantageous to substitute the time dependence of Ci(t) 
by time dependence of Ni(t); i.e., the concentration is substituted by the number of 
i objects in radical cluster, which is also in agreement with the fact that the number 
of SSB formed in DNA should be proportional to the number of radicals in cluster. 
Substituting Ci(t) by iVj(t) with the help of Eq. (16) the Eq.(18l may be rewritten 
as 



dNj 
dt 



= -N 



NjNj 
V 



VjV k 



(19) 



which is again the system of ordinary differential equations. 

The functions N(t) may be determined from the last equation (19) if the indivi 
dual Vi(t) functions may be regarded as given, which may be derived from Eq. ( 15 1 
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The volume Vi(t) may be then given as 



V t (t) = -Trff (t) 

where the average radius of the i species cluster equals 

1 



Ni(t) 



r Ci (r, t) Airr dr; 



(20) 



(21) 



and inserting according to Eq. ( 15 1 one obtains 



Vi{t) 



256 fDft 3 \ 



(22) 



From Eqs. (19) and (22 1 one obtains then the final system of ordinary differential 
equations 



dNi _ 3 
~dT ~ ~256 V Vi 3 



DPI) 



(23) 



Radical cluster and DSB formation 



Before passing to the analysis of experimental data we must distinguish between 
two kinds of species i = l,..,n involved in chemical reactions. Some of them, i.e., 
i = 1, .., s < n, are assumed to be radicals formed by the impact of ionizing particles, 
while the other are formed in the course of diffusion or are present permanently in 
the water solution. Thus the radical cluster is formed by the first kind of species 
that decrease quickly as the consequence of corresponding chemical reactions. 

It is evident that the actual origin of the cluster cannot be identified with t = 
of cluster volume evolution described by Eq. (15 1. It is necessary to introduce 
the time to that corresponds to initial values iVj(to), i = 1 5 --jS, forming a radical 
cluster immediately after the impact of ionizing radiation; the radicals reacting then 
with DNA molecules and forming individual SSB. However, according to our earlier 
results [10] an efficient cluster may be formed also in a certain distance from a DNA 
molecule; contact being realized later due to cluster diffusion and DNA molecule 
motion. Thus the probability of SSB formation by a radical i will equal 



Pi 



aiNi(t)dt, 



(24) 



where parameters cti (i = l,...,s) must be proportional to the reaction rates of 
individual species i with DNA molecules given in literature. As to the final time 
value it might be put t m — ► oo; however, practically it is sufficient to integrate to t m 
corresponding to Ni(t m ) < 1. 
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And for the probability p s of SSB formation in a DNA molecule it holds then 

Ps = /.Pi+Pdir, (25) 

i 

where parameter pdir represents a direct effect of the given ionizing radiation that 
might be additional to the indirect chemical mechanism. However, we are interested 
in radical clusters being able to form DSB in corresponding DNA molecules. The 
corresponding probability may be then given by 

Pd=P 2 s , (26) 

which may be correlated to experimentally established values of DSB numbers under 
different conditions (oxygen content in our case). 



Specification of the model according to data kind 



As already mentioned we shall analyze the data concerning DSB formation under 
different oxygen content. And we shall assume that the corresponding average radical 
cluster at time to consists of the following species: H', OH' 



The other species playing the role in final effect are then introduced in Table 
together with chemical reactions being considered in the right-hand side of Eq. d23 1 
the reaction rates (in corresponding units) established experimentally are also given 
Only a part of reactions introduced in Table [T] have been taken into account. 



Reactions 


Rate constants 

(nm a . N _1 . ns" 1 ) 


1. 


H' + H' 




H 2 


16.606 


2. 


eT aq + H' 




H 2 + OH- 


41.514 


3. 


f~ ~\- P~ 

^aq 1 ^aq 




H 2 + 20H- 


9.963 


4. 


e" + OH' 




OH~ + H 2 


49.817 


5. 


H' + OH' 




HO 


39.854 


6. 


OH' + OH' 




H 2 2 


6.642 


7. 


H a O + + e" 




H' + H 2 


38.193 


8. 


HO' 2 + H' 




H 2 2 


16.606 


9. 


Zaq + 2 




2 + H 2 


31.551 


10. 


HOI + OH' 




HO + 2 


16.606 


11. 


HO' 2 + HOI 




H 2 2 + 2 


0.003 


12. 


H' + 2 




HO' 2 


16.606 


13. 


2 + HO+ 




HO' 2 


49.817 



Table 3: Recombination reactions chosen from Table Q] 

The corresponding system of ordinary differential equations, which describe dy- 
namic of recombination and diffusion processes in the considered case may be then 
written in the form 



dN H 
dt 



' 256 



/(I 
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+ ieM-^- f W V**OMOJ + faO (27) 

256 V U ^j{D H30 D e f 

cLNqh = SN OH //tt\ / fc 5 jV g 2fc 6 7V og fc 10 jV o fc 4 7V e \ 
dt - - 256 VW^,/(D£) yp^y VCDS) " 

-iV OH (fc 14 [O a ] + fc 16 ) (28) 



3iVe 

' 256 



k 9 N 



k 4 N OH , 2fc 3 7V e 



H3O 



-fc 9 7V e [O a ] (29) 




2fcnJV 



fci 2 iV ff [Oa] - 



k 17 N 



(30) 



df 



256 Vvt 3 



k 7 N e 

7W) 



h 3 N r 



(31) 



dN r 



dt 



3N r 



256 



(?) 



ki 3 N H3 



Pko) 



k 9 N e [0 2 



(32) 



where ki, . . . , £;i3 are rate constant of reactions given in Table [3] and fci4, . . . , fci7 
correspond to additional reactions the rate constants of which were not available; 
they are taken as fully free parameters in the optimization procedure. On the other 
side the values of involved diffusion coefficients have been fixed as already mentioned 
(see Table |). 

In agreement with general opinion we will assume, that for SSB formation practi- 
cally the radicals H', OH',e~ q and HO' are responsible. The probability p s of SSB 
formation in DNA molecules depends then on the instant of encounter of a molecule 
with a given cluster. One can write for an average cluster 



PS = Ph +POH +Pe +PO+ Pi, 



(33) 
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where the individual probabilities arc given as averages over different t: 



PH = I a H N H (t)dt (34) 

to 



Poh = / aoiiNoH (t) dt (35) 

Jto 



Pe= f a e N e (t)dt, (36) 

Jto 



Po= / a N (t)dt (37) 

Jt 

and parameters an, ctoH, olo, ct e are held to be proportional to the reaction rates 
of individual radicals and of e~ q with DNA molecules (see, e.g., [5]). The value 
t m should correspond to cluster diffusion time (e. g. to time when less than one 



radical of the given type is present in a given cluster) . The parameter p\ in Eq. ( 33 1 
represents direct effect under the given irradiation conditions. 
The probability of DSB formation may be then given by 

Pd=Ps- (38) 



Analysis of experimental data 

As mentioned we have used already the preliminary version of the given mathemati- 
cal model in analyzing experimental data of Blok and Loman [3]. Even if their paper 
was published in 1973 the data represent still very important information about the 
formation of DSB (and also SSB) at different oxygen concentrations. The data were 
gained by irradiating the <1>X174-DNA in water solution by photons (~1,25 MeV) 
of Co-60 isotope; the applied dose was 5 Gy. The solution contained 25 fig.ml -1 
of DNA in 0,01 phosphate buffer at pH7. The main measurements concerned SSB 
numbers at different oxygen concentrations while the corresponding DSB numbers 
were established in some cases only to determine the ratio of DSB and SSB. We have 
established the corresponded dependence of DSB numbers on oxygen concentrations 
from the ratio of DSB and SSB numbers in all measured points (see Fig. 1). We 
have been interested in DSB formation mainly as they are responsible for biological 
effect of ionizing radiation. While the experimental data have contained the results 
for two different gas mixtures, in the following the model will be demonstrated being 
applied only to the mixture N2 — Oi- 

We have assumed that the DSB may be created in the case only when a greater 
radical cluster (formed during irradiation) meets a DNA molecule. And we have 
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Figure 1: DSB numbers (per DNA molecule) at different oxygen concentrations in 
N2 — O2 mixture; experimental data - individual points, theoretical dependence - 
full line. 



been looking for what the size of average cluster may be and which chemical pro- 
cesses may be running inside such a cluster before it diffuses. We have required for 
all chemical reactions being involved to possess chemical reaction rates correspond- 
ing approximately to values (or rather to their mutual ratios) being introduced in 
literature. And we have asked also for the ratios of initial numbers or radicals Nh , 
Noh 1 &aq t° correspond to frequencies of processes leading to their formation, i.e., 

N H - : Nqh : N e = 1:2:2. 

The initial values of other radicals have been put to be zero at time to as they may 
be formed with some delay after the original water radicals meet an O2 molecule. 
The agreement with literature data (with their mutual ratios) was asked also for re- 
action rates of different radicals with DNA molecules; i.e., for the ratios of quantities 
a H, ot e , ao ■ Only the rest of parameters have been left free to be determined 

in optimization. 

The actual values of individual parameters have been determined with the help 
of the MINUIT program (see [5])- The best fit to the given experimental data was 
obtained on the basis of the following values: For the initial numbers of radicals 
forming average cluster at to it has been obtained 

N H - = 18.11, Noh = 31.9, N e = 35.0. 

And for the reaction rates of different radicals with DNA (comp. [5] ) 

ot-H '■ oloh ■ a e : ao — 5.95 : 93.53 : 0.01 : 0.58. 
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The probability of a direct effect has been approximately 

px = 0.0003 , 

indicating to be practically negligible. The values of parameters ki, fci3 have been 
introduced in Table 3. And for the additional reaction rates it has been obtained 

fci 4 = 325.54, fci 5 = 0.098, k 16 = 0.107, k 17 = 0.004. 

The model dependence is shown in Fig. 1 together with the corresponding experi- 
mental data. 

Experimental curve exhibits, however, a quite flat dependence for higher values 
of oxygen concentration. It follows from the fact that the oxygen does not dissolve 
fully in water at all concentrations; the saturation depending on concentration value. 
It is then necessary to write for oxygen concentration [0 2 ] i n water solution (see Eqs. 
( 27p2| ) 

[0 2 ] = A {l-exp(-(ca ) T0 )} (39) 

where c is the oxygen concentration in gas mixture. Our fit has been then obtained 
when this dependence has been added with the following values of the involved 
parameters 

A = 0.0004, cr = 43.68, r G = 1.7; 

the value of parameter Ao representing maximal solubility in the given solution 
seems to be in good agreement with the value 0.0008 nm~ 3 for pure water. 

Conclusion 

The described model enables to study the influence of corresponding substances and 
of individual processes running during the chemical stage on the important damages 
of DNA molecules when irradiated by ionizing radiation. These processes may play 
important role especially for low-LET radiation (at not very high doses and dose 
rates) when individual DSB are formed mainly by track ends of secondary electrons. 
The approach might be easily adapted also to track ends of protons (or similarly) if 
the initial form of radical clusters is assumed to be cylindrical (instead of spherical) 
and the influence of diffusion is modified in corresponding way. 

In the presented paper we have demonstrated the influence of different oxygen 
content on DSB formation and the fact that the rather drastic experimental de- 
pendence in the region of very small oxygen content may be reproduced when all 
chemical reaction and diffusion rates are in agreement with ratios between corre- 
sponding values established in other experiments. Here we have tried to describe 
mainly the corresponding mathematical model and its solution approach without 
going to conclusions concerning proper radiobiological aspects. These goals will be 
followed in another paper. 

The model has opened, however, some other possibility of analyzing the influence 
of different contents of other present substances that may act also as radiomodifiers, 
enlarging or diminishing the radiobiological effect. When the reaction rates between 
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individual substances are known their effect at their divers concentrations and also 
at different oxygen concentrations might be established. A scries of corresponding 
experimental data that might be analyzed with the help of the presented model may 
be found in literature. 
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